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Scaling symmetries of the Euler-Lagrange equations are generally not variational symmetries of 
the action and do not lead to conservation laws. Nevertheless, by an extension of Noether's theo- 
rem, scaling symmetries lead to useful nonconservation laws, which still reduce the Euler-Lagrange 
equations to first order in terms of scale invariants. We illustrate scaling symmetry dynamically and 
' statically. Applied dynamically to systems of bodies interacting via central forces, the nonconserva- 

^ | tion law is Lagrange's identity, leading to generalized virial laws. Applied to self-gravitating spheres 
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I. SCALING SYMMETRY NOT GENERALLY A SYMMETRY OF THE ACTION 

Action principles dominate physical theories because they admit transformations among dynamical variables and 
exhibit common structural analogies across different systems. If these transformations are symmetries of the action, 
then by Noether's theorem, they give rise to conservations laws that reduce the number of degrees of freedom. 
This relationship of symmetries of the action (variational symmetries) to conservation laws is central to Lagrangian 
dynamics. Even if these transformations are not symmetries of the action, they nonetheless lead to a useful Noether's 
identity. 

Although equations of motion do not require Lagrangian expression, we apply Noether's identity to transformations 
that are not symmetries; in particular, to scaling symmetry, which is not generally a symmetry of the action, but only 
of the equations of motion (Section II). Variational symmetries and generalized symmetries both reduce the equations 
of motion to first order, but in different ways. 

• Variational symmetries imply conservation laws, first integrals of the equations of motion. 

• Scaling symmetry generally implies only a nonconservation law, which still reduces the equations of motion to 
first order in scaling invariants. 

Applied to dynamical systems of bodies interacting via inverse power-law potentials, these nonconservation laws 
are Lagrange's formulae, or generalized virial theorems (Section III). Applied to self-gravitating barotropic spheres 
in hydrostatic equilibrium (Section IV), the nonconservation law leads directly to a first-order equation for homology 
invariants and to well-known properties of polytropes and of homogeneous stellar cores (Section V). In this way, these 
nonconservation laws illuminate the physical consequences of scaling symmetry. 

The applications of continuous local symmetries of classical Lagrangians considered here are actually few. We 
do not consider applications to quantum field theories involving the symmetry of the vacuum state as well as 
the Lagrangian, which lead to important quantum anomalies and topological symmetries, generated by topological 
charges. 

II. NOETHER'S THEOREM EXTENDED TO SCALING SYMMETRIES 

A. Noether's Identity Implies Both Conservation and Nonconservation Laws 

We consider a one-dimensional discrete dynamical system described by the Lagrangian density C(t, qi, q\) and action 
S = J C(t,qi,q\)dt, where the dot designates the partial derivative d/dt with respect to the independent variable. 
Under the infinitesimal point transformation S(t, qi), 5qj(t, qi) generated by St ■ d/dt + 5qi ■ d/dqi, the partial derivative 
and the Lagrangian transform locally as 

d8qi dSt 

5C = [St ■ d/dt + Sq t ■ d/dq t + dqi ■ d/dq\]C = CSt + (dC/dq^Sq, + (d£/dq t ) - q~] , (1) 

where the total derivative d/dt = d/dt + q\ ■ d/dqi. The Einstein summation convention is assumed. 
The canonical momentum and Hamiltonian 

Pi(t,q t ,q t ) := dC/dq\ , h(t,qi,4i) := q^dC/ 'dq\) - C (2) 

have total derivatives 

= dC/dqi-T>i , = C + T>i- q\ , 

at dt 

l( Pl Sq l ) = (dC/dq l -V l )Sq l +p i -^ , = (£ + V . . g. )St _ h . &. , (3 ) 

in terms of the Euler-Lagrange variational derivative T>i := dC/dqi — d{dC/ dqi) / dt. The Noether charge 

G := C ■ St+pi ■ (Sqi - q\8t) = -hSt + pi8 qi (4) 
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has total derivative 



dG 
~dt 



SC + C 



d{8t) 
dt 



Di ■ (Sqi - q\8t) =8C — T>i ■ (Sq t - q\5t) , 



(5) 



where SC := SC + C- (d5t/dt) is the change in Lagrangian at a fixed point. Different Lagrangians, leading to the same 
equations of motion, define different Noether charges and nonconservation laws. 
The variation in action between fixed times and end points is 



SS 12 



dt SC 



dt 



dG_ c 
dt 



d{St) 
dt 



+ V i -{Sq i -q i St) = G(1)-G(2) + 





fdh 


dC\i 


c* 


Sqi-Vi+St-^—- 





(6) 



after integrating the term in d(St)/dt by parts. The action principle asserts that this variation vanishes for independent 
variations Sqi, St that vanish at the end points. It implies the Euler-Lagrange equations and dh/dt = —dC/dt. On- 
shell, where the equations of motion T>i = hold, 



SSu = J SCdt = G(l) - G(2) 



SC = 



dG 
~dt 



(7) 
(8) 



The last result is Noether 's equation, identifying the total derivative of the Noether charge with the change in La- 
grangian that generates it. By expressing the equations of motion in divergence- like form, it has an important 
consequence: a conservation law, if G is a variational symmetry; a nonconservation law otherwise. 

B. Transformations That Are Not Symmetries Still Lead to Useful Nonconservation Laws 



To the nonrelativistic central-force system 

h(r,p)=K + V(r):=p 2 /2m + V(r) , p 2 =p 2 r + l 2 /r 2 , 
apply the static radial dilatation (A) Sr = r and radial translation (B) Sr = r 1 := r/r: 

(A) : G := p • Sr = r • p: <5r = r , <5p = p , Sh = -2K + r(dV/dr) 

(B) : G := p • Sr — p r : Sr = ri , Sp = — Pt/r = r x 1/r 3 , Sh = —I 2 /mr 3 + dV/dr 
where 

p r = (r-p)r/r 2 , p f =p-p r = -rxl/r 2 , l:=rxp 



(9) 



(10) 



are the radial and transverse linear momenta, and the angular momentum, respectively Q. 

Because both these radial transformations are static, G = p ■ Sr, dG/dt = SC = —Sh, the two nonconservation laws 
are 

(A) : dG/dt = d(r ■ p)/dt = 2K - r(dV/dr) 

(B) : dG/dt = dp r /dt = l 2 /mr 3 - dV/dr . 

Except for circular orbits, neither of these transformations is a symmetry. Nonetheless, each of these nonconservation 
laws expresses important consequences of the equations of motion, in any central- force system @. 

(A) : Defining the virial p • r := A, the r • p nonconservation law is Lagrange's formula p = — W, which preceded 

Clausius by almost a century 0]. In the form dA/dt = 2K + r i ' Fj> the law still holds for a system of bodies, 
even if the central forces p,; = F, do not derive from a potential. 

(B) : The p r nonconservation law is the radial equation of motion. 

Both these nonconservation laws express the equations of motion and do not depend on scaling symmetry. But, if 
the potential is homogeneous in r, so that r(dV j ' dr)=—nV , the system is scaling symmetric. In any bounded ergodic 
system, the time averages (dA/dt) = (dp r /dt) vanish, so that 



4 



(A) : 2(K) = -n(V) 

(B) : (I 2 /m)(l/r 3 ) = —n{V/r) , 1^0 . 

For n = 1, (B) is useful for relativistic corrections to noncircular hydrogenic orbits. (A) is the usual virial law. 

In nondegenerate perfect gases, equipartition makes the internal gas kinetic energy K = | J PdV— J edm, where 
the internal energy density e = |(P/p) = %RT/n for a gas of molecular weight fi. The Coulombic virial theorem 
2(K) = —{V) then determines the averaged specific temperature (T//z) and leads to important applications in classical 
kinetic theory and in stellar structure. 

C. Variational Symmetries Imply Conservation Laws 

The most general and important applications of Noether's identity are to variational symmetries and to dynamic 
scaling symmetries of the equations of motion, which preserve the stationary action principle SS±2 — and reduce 
the equations of motion to first order in different ways. 

Variational symmetries preserve the action SSi2 — because SC — or dB/dt, the total derivative of some gauge 
term B(t,q). Noether's identity (d/dt)(G — B) = —T>i ■ (Sqi — q\5t) conserves G ~ B on-shell, when the equations of 
motion hold. This original version of Noether's theorem, identifying conservation laws with variational symmetries, 
has two familiar applications. 

Point symmetries lead to integration by quadratures: Any central-force system ([9]) is symmetric under time 
translations and spatial rotations, leading to conservation of energy E and angular momentum I: 

E := (r 2 +r 2 e 2 )/2 + V(r) , l:=mr 2 . (11) 
Since rO/f = p t jp r = v/ r 2 p 2 /l 2 — 1, quadrature leads to the first-order orbit and time equations 



8(r)=8 Q + J dr/r^2mr 2 [E -V(r)]/l 2 -1 , t = t + J rdr/ v / 2r 2 [E - V{r)\ - (l/m) 2 . (12) 

In the Newtonian case V(r) = —GM/r, the integrals reduce to elementary functions and the orbits are conic 
sections 

r(0)=p/[l-esin(0-0 o )] (13) 

of eccentricity e 2 =l + 2E(l/GmM) 2 , where p := (l/m) 2 /GM 0. 

Conservation laws including any gauge terms: A variational symmetry in which the Noether charge is not con- 
served obtains in the many-body system of particles with intcrparticle forces that depend only on the relative 
separations — q 3 and relative velocities q^ — cjj . This system admits the infinitesimal boost transformations 

5qi=5v-t , St = , SV = , SK = 5C = P ■ Sv , (14) 

where M,P,K are the total mass, momentum, and kinetic energy. The Noether charge G = (P • v)t is not 
conserved, but Noether's equation gives the conservation law (P — MR) • Sv = 0, or MR = P, for arbitrary 
infinitesimal Sv. Boosts change the total momentum P, but the center-of-mass moves with velocity R. This 
familiar center-of-mass theorem follows directly from boost symmetry, irrespective of the internal forces. It is 
paradigmatic for distinguishing between the effects of internal and external forces on many-body system. 

The converse of Noether's theorem is that conservation laws imply invariance of the Lagrangian up to a possible 
gauge term. For example, the conservation of the relativistic momentum and energy implies Lorentz invariance of the 
Lagrangian. 

D. Scaling Symmetry Implies a Special Nonconservation Law 



In a many-body system with individual coordinates r^, the dynamical scale transformation 

St = /3-t , Sri = Ti , 5{d/dt) = -(i-(d/dt) , STi = (l-/3)-Ti , S(r?/t) = (15) 
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is generated by the Noether charge 

G :=-/?• ht + A , (16) 

where A := X)P* ' r « * s the virial. Dynamical scaling is a symmetry of the equations of motion (but not of C), if 
the pairwise potential energies are inverse powers Vij ~ |r^ — r_,-|~ n of the interparticle distances, the potentials are 
homogeneous in their coordinates r(dVij/dr) = —nVij,SVij = —nVij and (3 = 1 + n/2, so that all distances scale as 

n ~ tW = t»/( 2 +«) 0. 

Scaling symmetry makes the Lagrangian a homogeneous function of its arguments, a scalar density of some weight 
— 2cj, so that SC = —2CoC, SS12 = oS\2, where a := 1 — 2ui. Noether's identity then implies the special on-shell scaling 
nonconservation law 

An 

^ = (l-2u)C = aC , (17) 

which reduces to a conservation law only asymptotically, wherever aC is small. This asymptotic conservation law 
then allows approximate integration of the equations of motion, in certain limits. 

In the next section, we consider energy-conserving mechanical systems £=£(r,-, r,), for which the dynamic scaling 
nonconservation law is a generalized virial law. In Section IV, we consider the spherical hydrostatics of barotropic 
fluids, for which the Lagrangian C(r, H, H') depends explicitly on the independent variable r. Instead of a first 
integral, both these examples illustrate a first-order differential equation among homology invariants 0, @ , linearly 
relating the "energy function" h = K + V to C = K — V, or the "kinetic" term K to the "potential" term V. 



III. DYNAMICAL NONCONSERVATION LAWS FOLLOWING FROM SCALING SYMMETRY 

A. Mechanical Nonconservation Laws Are Generalized Lagrange Identities 

Consider a nonrelativistic system of particles with coordinates r,, momenta Pi, interacting by pairwise static 
potential energies Vij . The dynamical scale transformation (|15p generates the infinitesimal changes 

5K = 2(1- /3)K , 5V = -(l-/9)r- W , SA = (1-/3)A . (18) 

If the pairwise potential energies are inverse powers Vij ~ |r^ — rj |~™ of the interparticle distances, the potentials are 
homogeneous in their coordinates rij(dVij / 'rfry) = —nVij,5Vij = —nVij. Provided n = 2(/3 — 1), j3 = 1 + n/2, the 
Lagrangian density scales as 

SC = -nC , 5C = 6C + f3C = (1 - n/2)C . (19) 

Because energy is conserved, dG/dt = —(3h + A, so that the scaling symmetry nonconservation law 

A = (l + n/2)h+(l-n/2)C = 2K + nV (20) 

relates the nonrelativistic kinetic energy K and power-law potential V to the time derivative of the virial A. 

For periodic or long-time averages in bounded ergodic systems, we have (A) = 0, and the virial theorem 2(K) 
= — n(V). Table I tabulates these generalized virial theorems and period-amplitude relations for orbits in the five 
important inverse-power-law potentials n = —2,-1,0,1,2. Only for n = 2 does dynamical scaling reduce to a 
variational symmetry, so that the Noether charge G = —2(K + V)t + A is conserved. For potentials more singular 
than 1/r 2 , there are no bound states. 



TABLE I: Period- Amplitude Relations and Virial Theorems for Inverse Power-Law Potentials V ~ 1/r 



n 


System 


Period-amplitude relation t ~ r 1+n / 2 


Virial theorem 


-2 


isotropic harmonic oscillator 


period independent of amplitude 


(K) = 


(V) 


-1 


uniform gravitational field 


falling from rest, e.g., z = gt 2 /2 


(K) = 


(V)/2 





free particles 


constant velocity r ~ t 


(K) = 





1 


Newtonian potential 


Kepler's Third Law t 2 ~ r 3 


(K) = 


-<V>/2 


2 


inverse cube force 


t ~ r 2 


(K) = 


~(V) 
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B. Scaling Nonconservation Law in Classical Electrodynamics 

Noether's identity applies to continuous Lagrangian systems (fields) as well as discrete systems. In this case, r, t 
are independent variables. If f , G = E x B/47rc, T, U are respectively the electromagnetic force density, momentum 
density, momentum flux tensor, and energy density, then momentum balance reads 

dG/dt + V -T+i = . (21) 

From this follows an electromagnetic analogue of the mechanical Lagrange's identity: 

<3(r-G)/5t + V-(T-r)-[/ + r-f = . (22) 

When time-averaged, this becomes an electromagnetic virial theorem Q. 



C. Scaling Nonconservation Law in Classical Conformal Field Theory 

In any relativistic field theory, space-time scaling (dilatation) symmetry leads to the familiar nonconservation law 



dx^ 



= ©£ , (23) 



where d/dx^ is the four-dimensional divergence, G M is the dilatation current, and 0£J is the trace of the energy- 
momentum tensor [H 0j @j- The dilatation charge is conserved when this trace vanishes, implying a conformal 
symmetry. 

The most familiar example of conformal symmetry is Laplace's equation in n spatial dimensions. In two dimensions, 
conformal symmetry implies the Cauchy-Riemann equations, so that any analytic function is a solution of Laplace's 
equation. In higher dimensions, conformal symmetry implies the conservation laws associated with translations, 
rotations, dilatations, and spatial inversions. The pure electromagnetic field is conformally invariant. 

These dynamical systems illustrate how Noether's identity leads to useful and often familiar nonconservation laws, 
even when scaling symmetry is broken. The remainder of this paper considers the hydrostatic equilibrium of gaseous 
spheres, where the independent variable is the radial coordinate r and the variational principle is that of minimum 
energy in place of least action. 



IV. SCALE-INVARIANT BAROTROPIC SPHERES 



The structure of luminous stars depends on the coupling between hydrostatic and thermal structure through an 
equation of state P = P(p, T), which generally depends on the local temperature and chemical composition. Ignoring 
chemical evolution, the matter entropy is locally conserved, so that in the steady state, stars are in both local thermal 
and chemical equilibrium. To treat the hydrostatic equilibrium independently of heat flow, we consider only stars in 
which the thermal structure is specified independently, so that the local equation of state is barotropic P — P(p) and 
dP/p = dH in terms of the specific enthalpy H (r) = E + Pj p. This restriction to barotropic stars makes the density 
p(r), specific internal energy E(r), specific enthalpy H(r) = E + P/p, and thermal gradient V(r) := dlogT/dlogP 
implicit functions of the gravitational potential V(r). The assumption of the independence of the thermal structure 
is justified as a good approximation if, as is usually the case, the thermal (Kelvin-Helmholtz) relaxation time of the 
whole star is much longer than its hydrostatic equilibration time. 

The hydrostatic structure of barotropes depends only on two first-order equations, mass continuity and of pressure 
equilibrium, 

dm/dr = Airr 2 p , -dP/dr = Gpm/r 2 , (24) 

or 

= -^°P > (r 2 H'Y + 47TGr 2 p(H)=0 , (25) 

r z dr V p dr / 

where ' := d/dr. In terms of the specific gravitational force dV/dr = g := Gm/r 2 , the equation of hydrostatic 
equilibrium reads 

d(H + V)/dr = . (26) 
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FIG. 1: Effect of the outer boundary on the Emden polytrope density profiles. For the softest equation of state n — 5, the 
stellar radius is infinite, u>5 = (5/3)(3 — it) everywhere, and p = p% p 3//6 . For stiffer equations of state n < 5, w n ~ wb inside 
the core u < 2, but increases as the finite radius is approached u — > 0. 



These structural equations are the Euler-Lagrange equations of the Lagrangian 

C(r,H,H') =4Trr 2 [~(H') 2 /8TTG + P(H)} , (27) 

derived from a minimal energy variational principles in Appendix A. However, the following consequences of assuming 
scale invariance do not depend explicitly on a Lagrangian formulation or on Noether's identity. 



A. In a Simple Ideal Gas, Scale Invariance Requires a Constant Entropy Gradient 



Polytropes are barotropic spheres with constant polytropic exponent 1 + 1/n :— d log P/d log p and polytropic index 
n = d log p/d log H. The pressure, specific energy, specific enthalpy, enthalpy gradient and central pressure at any 
point are 

P/P c = (p/p c ) 1+1 / n , E = n(P/p) , H=(n + l)(P/p) , dlogH/dlogP = l/(n + 1) . (28) 
Polytropic mechanical structure does not fix the thermal structure, which depends on the heat transport mechanism. 
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In a simple ideal gas, the equation of state, specific internal energy, specific enthalpy and adiabatic exponent are 

P/ P =-T , E = C V T , H = C P T , {dlogP/dlogps , (29) 
A* 

where 1Z is the universal gas constant, p is the molecular weight, and Cp, Cy are the specific heats at constant pressure 
and density. Following the law of energy conservation (Appendix B), the specific entropy and thermal gradient of a 
simple ideal gas are 

dS = C v d log P-Cpd log p , S = C v log (P/p^) , V = dlogH/dlogP . (30) 

Bound in a polytropic structure of index n, an ideal gas has constant thermal gradient, gravithermal specific heat, 
and entropy-pressure gradient: 

V = l/(n + l) , C* =C P (1-W ad /W) , dS/d\ogP = C P (V- V ad ) . (31) 

The radial entropy gradient 

dS/d log r = C P ( V ad - V) • v n , (32) 

is proportional to the homology invariant v n and positive (zero) when the thermal gradient is subadiabatic (zero). 
(See Section IV. B for more about v n .) In convective equilibrium, any polytrope has constant entropy 5*. In radiative 
equilibrium, a simple ideal gas polytrope has constant temperature and entropy gradient dS/ d log P. 

The thermal gradient is nearly constant and the hydrostatic structure nearly polytropic in zero-temperature (de- 
generate) stars and in chemically homogeneous stars starting out on the hydrogen-burning, zero-age Main Sequence 
(ZAMS): 

White dwarfs and neutron stars: nonrclativistic and extreme relativistic degenerate stars; polytropes of index 
n=3/2 and 3, respectively. 

ZAMS stars in convective equilibrium: with vanishing gravithermal specific heat C* and uniform entropy den- 
sity. These are n=3/2 polytropes. 

ZAMS stars in radiative equilibrium: With uniform energy generation and Kramers opacity, stable polytropes 
of n > 3/2. At zero age, our Sun was a chemically homogeneous star of mean molecular weight p = 0.61, well- 
approximated by the Eddington standard model (n=3) throughout its radiative zone, which contained 99.4% of 
its mass. Because energy generation was centrally concentrated, our ZAMS Sun was even better approximated 
by a slightly less standard n=2.796 polytrope @. 

Even better nonpolytropic fits, to the M—R relation R ~ observed in young ZAMS stars are obtained by including 
nonuniform energy transport and corrections to Kramers opacity: radiative transport in the pp-burning lower 
main sequence 0.11 < M/Mq < 1.2 gives £ = 0.57; convective transport in the CNO-burning upper main 
sequence 2 < M/M Q < 20 gives £ = 0.8 [13, EH- 

Chemically inhomogeneous stars and the photospheres of luminous stars cannot be polytropic. Because our present 
Sun is chemically evolved and has a convective envelope, it is far from being polytropic: its polytropic fit, with index 
n = 3.26, is poor @. 

B. Polytropic Structure Implies a First-Order Equation in Scaling Invariants 

Following Chandrasekhar [ijj], we define homology variables 

u := d log m/d log r = 3p/p , v :— ~d\og (P/p)/d\ogr , w := —dlog p/d\ogr — n(r) ■ v , (33) 

where p — 3m/4irr 3 is the average mass density interior to radius r and n(r) := dlog p / dlog (P/p). The central 
boundary condition is 

u(0) = 3 , u(0) = , (dv/du) = -5/3n . (34) 
The mass continuity and hydrostatic equilibrium equations (|24j) become 

dlog u/d log r = 3 — u — n(r)v , dlogu/dlogr = u — 1 + v — dlog [1 + n(r)]/d log r , (35) 



0.5 1 1.5 2 2.5 



: d logm/ d log t = r iipj~p 



FIG. 2: Dilution of polytrope density profiles in the envelope as the boundary is approached (u — ► 0). All solutions approach 
the same density structure w n (z) — s- ws = (5/3) (3 — u) at the center (u — > 3), but differ outside the core. As the boundary is 
approached, the mass m — > M, and the density scale height r/TO n — >• 0. 

are autonomous only when the index n(r) is constant. 

In polytropes, the constant index n and gradient V = l/(n + 1) imply 

du/dlogr = u(3 — u — nv n ) , dv n /dlogr = v n (u — 1 + v n ) , (36) 

which are both autonomous and can be written as the characteristic equations 

d log v n d log u 



U — 1 + V n 3 — U — nv n 



= d\ogr (37) 



for the phase variables u, v n . Deferring the second equality for r = r[v n (u)] to the next section, we now solve the 
first-order Abel equation 

dv n v n (u-l + v n ) 

~r~ = ~n v ' ^> 

au u{A — u — nv n ) 

subject to the central boundary condition u(0) = 3,v n (0) = for regular (Emden) polytropes v n (u) [l7l ]. 
The infinitesimal scale transformation 

5r = r , 5H = -QH , 8H' = -(I + oj)H' (39) 
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is generated by the Noether charge 
whose radial derivative 



R ,2 H n+1 

r 1 

V 2 n + 1 



■l>66' 



r z , (40) 



dG n /dr = 6£ + £+V r -Q- d{Hr) jdr , (41) 
obeys the scaling nonconservation law 

dG n /dr = aC , a:=l-2u> , (42) 

wherever the Euler-Lagrange equation (25) is satisfied. The structural equation is scale-invariant if and only if 

■mo = 2 + o) , 2(1 + u>) = (n + l)u> , u = (n - 5)/(n - 1) , (43) 

so that the Lagrangian (27) is homogeneous of degree — 2w and <5£ = —2ujC. The scaling nonconservation law (jl!?)) 
then connects the gravitational and internal energy densities, just as the point-mechanics Lagrange identity connected 
the potential and kinetic energies. 

For polytropes, we introduce dimensionless units 

e:=r/a , 9 n :=H/H c = (p/p c ) 1 / n , (44) 



and the dimensional constant 



a 2 := { -^-K P l/ n - 1 = (n + l)/4nG-(P c /p 2 c ) , (45) 



where p c is the central density and 

P c /p c :=Kp 1 J n = 1 + n , JI :=(n + l)P c /p c = (n + l)li> 1 / n . (46) 
The included mass, mass density, average mass density, and gravitational acceleration are 

m = 4w* 3 • (-CX) , P = Pc-K , p = p c -(-39'JO , <7 = 4 W * 2 (-0 ■ (47) 

The Euler-Lagrange equation (f2"5"j). combining mass continuity and hydrostatic equilibrium, takes the dimensionless 
Lane-Emden form 

!( ? *) + «_0 . (48) 
In terms of 0, 9' , the homology variables [l2[ 

u:=dtogm/dlogr = 3p(r)/p = -#£/0J, > v n := -d\og(P/p)/dlogr = -£9'J9 n , p/p c = -3^/^ 2 , 

u/v n =q;+ i /o2 , ««„ = a?r x = r^r- 1 , w»:=^ +1 (-o = («<)* /a • ( 4 9) 

Hereafter ' := d/d£, and polytropes of different index are distinguished by the subscript n attached to different 
homology variables. 

Extracting the dimensional constant C := —H 2 /G := —[(n + 1)K pV""] 2 /G and suppressing the subscript n on 9, 
the Lagrangian, Hamiltonian, and Noether charge are 

£/C = £ 2 [9' 2 /2-9 n+1 /(n + l)} , 
n/C = e[9' 2 /2 + 9 n+1 /(n + l)] , 

G n /c = e d^ + — T ) +s>ee f . (50) 

L V 2 n+1/ J 

The scaling nonconservation law ((32) 

■ K^)+HH<?-^) 



is equivalent to the Lane-Emden equation (l48l) . It describes the evolving ratio between local internal and (negative) 
gravitational energy densities 

9 n+1 /(n+l) _ 2 u 

9 r2 /2 n + lv n ' 1 j 

as the local energy density changes from entirely internal at the center, to entirely gravitational at the stellar surface. 



11 



C. Scaling Fixes the Mass-Radius Relation Characterizing Different Polytropes 

We consider only Emden functions, which are regular at the origin and normalized to 9„(0) = 1, 9' n (0) = 0. Their 
first zeros 0(£i n ) = determine the stellar radius R = a^\ n . Each Emden function of index n is characterized 
equivalently by its dimensionless outer radius £ ln , its outer boundary value u> n of the homology invariant oj n , or its 
density ratio p cn /Pm where p n = SM/AirR 3 is the mean density. All three are tabulated in the third to fifth columns 
of Table II, for eight values of the polytropic index n. Scaling relates the mass and radius, according to the M-R 
relation M ~ ow„i?' Tl_3 - ) ^ Tl_1 - ) in the last column. 

We define the inner core radius £i cn implicitly by w(£icn) = 2, the radius where the acceleration Gm/r 2 reaches a 
maximum and the gravitational energy density overtakes the internal energy density. The sixth and seventh columns 
in Table II list dimensionless values for this core radius ri cn /R — £i crl and included mass mi cn /M, shown by red dots 
in Figures 3, 4. According to equation (|52p. the internal energy dominates in the core; while in the envelope, the 
gravitational energy dominates. 

For homology variables, we prefer a new independent variable z := 3 — u = —d\ogp n /d\ogr and a new dependent 
variable w n :— nv n :— —d\ogp/d\ogr. In term of these variables, the hydrostatic equilibrium and of mass continuity 
characteristic equations (I3T)) are 

dz/(3 — z)(w n — z) = d\ogw n /(2 — z + w n /n) — dlogr . (53) 

The first equality is the first-order Abel equation for the invariant w n (z), which we solve for the central boundary 
condition w n — > 5z/3 when z — s> 0. w n (z) and the differences (3/5)[w n (z) — w^(z)] = (3/5)w n (z) — z = dlogp/ps ~ G n 
are plotted in Figures 2 and 1, respectively, for polytropic indices n = 1, 2, 3, 4, 5. 

For n < 5, the stellar boundary lies at finite radius. The Noether charge is nearly conserved at G n ~ in the inner 
core, but grows rapidly as the boundary radius is approached (Figure 1). In the envelope, the growing Noether charge 
measures how close the boundary is. The finite radius R determines the stellar scale, even though the polytropic form 
is locally scale invariant . 

For incompressible matter (n — 0), there is no core concentration: the mass is uniformly distributed, and the entire 
star is core. But as the equation of state softens as n increases toward 5, the gradient decreases, the core concentrates, 
the inner core radius shrinks, and the envelope outside the core grows: r- lcn /R — > 0, m,\ cn jM — »~ 0.19. For the softest 
equations of state n < 5, the stellar radius £i„ w 3(n + l)/(5 — n), the inner core radius shrinks £; cn « ^/l0/3n, their 
ratio r icn /R = £icn/£in « 0.045(5 - n), m icn /M w 0.20, and uj n w y/3/£i n . 

For n = 5, the core becomes infinitely concentrated, shrinking to zero, and the star is all envelope. The n = 5 
regular solution 

M£) = (l + e 2 /3)- 1/2 , p = p c (l + e/3y 5/2 , m = Me/(3 + e) 3/2 , « 6 = (3-«)/3 (54) 

has infinite stellar radius R as shown at the bottom of Table II. In this case, not only is the differential form scale 
invariant, but also the action and stellar structure. The n — 5 polytrope is globally scale invariant, and the Noether 
charge G5 ~ (v^ + u/3 — 1) = 0. 
For < n < 5, 

w n (z) = I" dz w n % - Z + Wn,n \ « (5/J)[l - (1 - z/3) J ] := w nPic (z) , J := (9n — 10)/(7— n) , (55) 
Jo \6 — z)[w n — z) 



TABLE II: Scaling Exponents, Core Parameters, Surface Parameters, and Mass-Radius Relations for Polytropic Gas Spheres 
of Increasing Core Concentration 



n 






pen j Pn 




Hen /R 


m ic n/M 


M ~ R 1 -"™ Properties 





-2 


2.449 


1 


0.333 


1 


1 


M ~ R 3 ; incompressible matter, all core 


1 


±00 


3.142 


3.290 




0.66 


0.60 


R independent of M 


1.5 


4 


3.654 


5.991 


132.4 


0.55 


0.51 


M ~ R~ S ; nonrelativistic degenerate 


2 


2 


4.353 


11.403 


10.50 


0.41 


0.41 




3 


1 


6.897 


54.183 


2.018 


0.24 


0.31 


M independent of R; Eddington standard model 


4 


2/3 


14.972 


622.408 


0.729 


0.13 


0.24 




4.5 


4/7 


31.836 


6189.47 


0.394 


0.08 


0.22 




5 


1/2 


00 


00 








0.19 


maximally concentrated; entirely envelope 
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FIG. 3: Normalized density profiles as a function of fractional included mass m/M, for polytropes of finite mass M and 
compressibility increasing with n. The red dots mark the cores. For incompressible matter (n = 0), the polytrope is all core. 
As the matter softens (n increases), an envelope grows to ultimately encompass just over 80% of the mass. For any n > 1, the 
density at the inner core radius stays in the narrow range 0.37 < p(ri cn )/p c < 0.42. 



TABLE III: Taylor Series and Picard Approximations 8 n pi c to Emden Functions 0. 



n 


Emden Function and Taylor Series 


N := 5/(3n - 5) 


Picard Approximation 0„pj c := (1 + £ 2 /6iV) N 







-1 


l-C 2 /6 


1 


sin£/£ = 1 - C 2 /6 + £, A /12Q - £ 6 /5040 + . . . 


-5/2 


(1 - C 2 /15) 5/2 = 1 - C 2 /6 + C 4 /120 - £ 6 /10800 + ■ ■ • 


n 


1 - £ 2 /6 + < 4 / 120 - n(8n - 5)/15120£ 6 + • • • 


5/(3n - 5) 


(1 + f/6N)~ N = 1 - £ 2 /6 + n£ 4 /120 - n(6n - 5)£ 6 /10800 + 


5 


(1 + C 2 /3)" 1/2 


1/2 


(i+em~ i/2 



is well-approximated by the Picard approximation obtained by inserting the core values w n (z) « (5/3)z inside the 
integrals. Indeed, this Picard approximation is everywhere exact for n = 0, 5. For < rt < 5, it breaks down only in 
the outer envelope, where w n diverges as w n —> n[ow" _1 /ii] ly/n , and w„ must be calculated from the exact asymptotic 
value of u> n given in Table II. 

Integrating over z, the density profile and Emden functions are Q 

Pn (z)/ Pcn = K = exp{- f dzW " (z) } « (1 - */3) 5 / 2 (56) 

{-/' F fr"hl , }^(l-^/3) 5/2 ":^^Pic , (57) 
I- Jo n[w n {z) — z\(3 — z) ) 



exp 



where again the Picard approximations are obtained by inserting the core relations w n (z) = nv n (z) ss (5/3)2 under 
the integrals. 
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r/R 



FIG. 4: Normalized density profile as function of fractional radius r/R. The density is constant for incompressible matter 
(n = 0) which is all core, but is concentrated at the origin for an unbounded polytrope (n = 5, R = oo) which is all envelope. 
For any n > 1, the density at the inner core radius stays in the narrow range 0.37 < p(ri cn )/p c < 0.42. 

V. CLOSED FORM APPROXIMATION TO EMDEN FUNCTIONS 

With the solutions w n (z) to the first-order equation, we now use the second equation (I3T1) 

dm/m :— u ■ dr/r = dz/\w n (z) — z] (58) 

to obtain 

,„ W /m = W3) v* . exp {^ d2 {_J_ _ £}} a (z/3) m (59) 

r(z),R - (,/»)* ■ exp 1»{ (3 _, ) ,J t( , ) - ^}} - M"7(3 - .) (60) 

for the mass and radial distributions. The integration constants i?, M, p c express the scale dependence of the 
polytrope. 

Using the radial distribution (58) to eliminate z(£), the Picard approximations 

0nPic(O = (1 + Z 2 /6N)- N , JV:=5/(3n-5) (61) 

to the Emden functions are obtained and tabulated in the last column of Table III. For n — and 5 polytropes, this 
closed form is exact. For intermediate polytropic indices < n < 5, the Picard approximation breaks down near the 
outer boundary, but remains a good approximation over most of the polytrope's bulk. 

Indeed, the Picard approximation is far better than any truncation of the Taylor series expansion of 8 n , whose 
radius of convergence is £ ~ 2. For the worst case, the Eddington standard model (n = 3), the Picard approximation 
$3Pic(£) to the exact Emden function and its tenth-order polynomial approximation: 

03(f) ~ 1 - £ 2 /6 + £ 4 /40 - (19/5040)<f + (619/1088640)£ 8 - (2743/39916800)£ 10 (62) 

are shown in Figure 5. Because this Picard departs from the Taylor series expansion already in sixth order 

03Pic(O = (l + 2f 2 /15r 5/4 = # 3 + e73528+--- , (63) 
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FIG. 5: The exact Eddington standard model Emden function 9a(£,), its Taylor series truncated at tenth-order, and its Picard 
approximation. Even in this worst case, the Picard approximation works well outside the core radius at £j C 3 = £i3(n C 3/-R) = 
1.65, but breaks down near the boundary. For nonstandard polytropes (n 7^ 3), the Picard approximation is even better and 
becomes exact everywhere as n — > or 5. 

it remains 90% accurate out to £ ~ 3.5, more than twice the core radius and more than half-way out to the stellar 
boundary at £13 = 6.897. Except for their very outer envelopes, which contain little mass and are never polytropic, 
the Picard approximations in white dwarf and ZAMS stars should be even better than for this n = 3 polytrope. 

VI. CONCLUSIONS 

We have generalized Noether's theorem connecting variational symmetries to conservation laws to generalized 
symmetries of the Euler-Lagrange equations. Although these lead only to nonconservation laws, they still reduce the 
Euler-Lagrange equations to first order, plus a quadrature. For scaling symmetries, the nonconservation law takes a 
special form linearly connecting the "kinetic" and "potential" parts of the Lagrangian. In special cases, a symmetry 
of the Euler-Lagrange equations and the nonconservation law can reduce to a conservation law and symmetry of the 
Lagrangian, the action, and possibly the solution. 

For nonrelativistic systems with inverse power law potentials, the scaling nonconservation law is a Lagrange's 
identity, leading to generalized virial theorems. For spherical hydrostatic systems obeying barotropic equations of 
state, the scaling nonconservation law leads to an analogous linear relation between the local gravitational and internal 
energies. From this nonconservation law, we derive all the properties of polytropes. Quadrature then leads to the 
regular (Emden) functions and their Picard approximations, which are useful wherever stars are approximately or 
exactly polytropic. 
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Appendix A: LAGRANGIAN FORMULATION OF BAROTROPIC HYDROSTATICS 

Stellar structure generally depends on coupled equations for pressure equilibrium and heat transport. Only if the 
heat transport leads to a local barotropic relation P = P(p) can the hydrostatic equations be considered independently. 
In such barotropes, the mechanical structure is fixed without reference to the thermal structure 

1. Mass Continuity and Hydrostatic Equilibrium 

We consider a self-gravitating isolated system in local thermodynamic equilibrium, a barotrope held at zero external 
pressure. The thermodynamic potential energy or work need to adiabatically extract unit mass is the specific enthalpy 
H(p) = E + P/p. Barotropic energy conservation, dH := dP/p, makes the specific enthalpy a more natural state 
variable than the specific internal energy E, pressure P, or density p. The equation of hydrostatic equilibrium 
—dP/dr = Gmp/r 2 := pg is then 

- dH/dr = dV/dr = g = Gm(r) /r 2 , (Al) 

describing how this local specific enthalpy or extraction energy H{r) depends on the local gravitational potential 
V(r). Integrating, we have the energy conservation equation 

H(r) + V(r) = -GM/R , r < R , (A2) 

where the zeros of the gravitational potential and specific enthalpy have been chosen at infinity and at the spherical 
surface, respectively. 

Because the gravitational potential obeys Poisson's equation 



1 d / dV 
r 2 d 

the specific enthalpy obeys the second-order equation 



\7 2 H + 4irGp(H) = . (A4) 



Implementing the equation of hydrostatic equilibrium requires a local entropic relation P(p), or P(H), p(H), which 
is determined by the thermal stratification of the static matter distribution in local thermodynamic equilibrium, and 
by a central boundary, or regularity, condition (dP/dr)o = = (dp/dr) . Near the origin, 

p^p c (l-Br 2 ) , m {r)^^^(l-*-Br 2 ) , 
p(r) := 3m(r)/47rr 3 -> p c (l - -Br 2 ) = p 2 c /5 p 3/5 , dw/du ->• -5/3 , (A5) 
in terms of the homology variables w := — dlogp/dlogr, u := dlogm/dlogr for the mass density and included mass. 

2. A Constrained Minimum Energy Principle for Hydrostatic Equilibrium 

In a static, self-gravitating sphere of mass M and radius R, the Gibbs free energy 

W :=E-TS + PV = n + U , (A6) 
in terms of the gravitational and internal energies 

fl = - (Gm/r)dm , U=- PdV , (A7) 
Jo Jo 

where p, E, and —Gm(r)/r are the mass density, specific internal energy, and gravitational potential, respectively. In 
the Eulerian description, the radial coordinate is r, the enclosed volume is V = 47rr 3 /3, and the enclosed mass m(r) is 
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constrained by mass continuity dm{r) — pdV. The Gibbs free energy is the work available in adiabatically expanding 
the sphere at fixed external pressure. If 



w = 



pit pit 

I £{r,m,m')dr = 4irr 2 [Gmp/r + P(p)]dr , (A8) 

o Jo 

abbreviating ' := d/dr, then the Lagrangian C is the Gibbs free energy per radial shell dr. 

The constrained minimum energy variational principle fill [l3| for hydrostatic equilibrium is that the Gibbs free 
energy be stationary (SW = 0) under adiabatic deformations in specific volume SV P = d(Aitr 2 8r) / dm that vanish 
on the boundaries and satisfy the mass continuity constraint m' — Airr 2 p. This minimum energy principle has the 
equation of hydrostatic equilibrium 

XV := Gm/r 2 + H' = Gm/r 2 + dP/pdr = , (A9) 

as its Euler-Lagrange equation, with mass continuity as a constraint. This equation is scale invariant if the specific 
enthalpy H scales as m! . 

3. An Unconstrained Variational Principle 

Using Poisson's equation to incorporate the mass continuity constraint, the gravitational energy is 

fi = - / (V' 2 /2)4irr 2 dr , (A10) 



so that the second-order Lagrangian (used in Section IV) 

C(r, H, H') = 4irr 2 [-H' 2 /87rG + P{H)] (All) 

is unconstrained and has Euler-Lagrange equation (|A4j) . The canonical momentum and Hamiltonian are 

p:=dC/dH' = -r 2 H'/G = -m , U{r, H,p) = -Gp 2 /2r 2 - Anr 2 P{H) , (A12) 

and the canonical equations are 

dH/dp = H' = -Gp/r 2 , dU/dH = -p' = m' = Anr 2 p . (A13) 

Spherical geometry makes the system nonautonomous, so that dH/dr = —dC/dr = —2C/r vanishes only at large r, 
with vanishing sphericity. 

Appendix B: STELLAR THERMODYNAMICS AND CONVECTIVE STABILITY 

The structure of luminous stars depends upon the coupling between hydrostatic and thermal structure through 
an equation of state P = P(p,T, p), which generally depends on the local temperature and chemical composition. 
But, ignoring evolution, the matter entropy is locally conserved, so that steady-state stars are in both local thermal 
equilibrium and mechanical equilibrium. In a fluid held in pressure equilibrium at constant external temperature, the 
specific Gibbs free energy H — TS = —V(r) is a minimum. Under hydrostatic equilibrium, the density p(r), specific 
internal energy E(r), specific enthalpy H{r) = E + Pf p, specific entropy and thermal gradient V(r) := d log T/d log P 
depend implicitly on the gravitational potential V(r). 

In the first law of thermodynamics 

TdS = dQ = C v dE + Pd{l/p) , (Bl) 

6 1 



E, p can be written as functions of temperature and pressure. Clever use of thermodynamic identities then leads 

to Ed, EH 



TdS = C*dT , dS = C P (V- V ad )dlogP , (B2) 



where the gravithermal specific heat C* := dS/d\ogT = Cp(l — V a£ //V) depends on the specific heat Cp and on 
the adiabatic gradient V a( 2 := (<91ogT/<91ogP)s. This expression of the first law of thermodynamics relates the local 
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thermal gradient V(r) := dlogT/dlogP to the gradient of the specific entropy S(r), which derives ultimately from 
the heat transport and generally varies in stars that are not in convective equilibrium. 

According to Schwarzschild's minimal entropy production criterion, convective stability requires dS/d\ogP < 0, so 
that the specific entropy is constant in convective equilibrium and increases outward in radiative equilibrium. This 
makes barotropic stars of mass M extremal in two respects: the central pressure is minimal in barotropic stars of a 
given radius R; the central pressure and temperature are maximal in barotropic stars of given central density. Because 
stellar evolution is driven by developments in the core, these bounds drive stars toward uniform entropy in late stages 
of evolution 
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